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Abstract. We study the random graph obtained by random deletion 
of vertices or edges from a random graph with given vertex degrees. 
A simple trick of exploding vertices instead of deleting them, enables 
us to derive results from known results for random graphs with given 
I*^ , vertex degrees. This is used to study existence of giant component and 

existence of fc-core. As a variation of the latter, we study also bootstrap 
percolation in random regular graphs. 

We obtain both simple new proofs of known results and new results. 
An interesting feature is that for some degree sequences, there are several 
or even infinitely many phase transitions for the fc-core. 



1. Introduction 



One popular and important type of random graph is given by the uni- 
K^ I formly distributed random graph with a given degree sequence, defined as 

^^ ■ follows. Let n G N and let d = (di)" be a sequence of non-negative integers. 

VQ '. We let G{n, d) be a random graph with degree sequence d, uniformly chosen 

among all possibilities (tacitly assuming that there is any such graph at all; 

t:;;}^ . in particular, ^^ di has to be even). 

CP I It is well-known that it is often simpler to study the corresponding random 

multigraph G*{n, d) with given degree sequence d = (dj)", defined for every 
sequence d with Yli ^i even by the configuration model (see e.g. Bollobas 
[3|): take a set of di half- edges for each vertex i, and combine the half-edges 

S^ . into pairs by a uniformly random matching of the set of all half-edges (this 

H \ pairing is called a configuration)] each pair of half-edges is then joined to 

form an edge of G*{n, d). 

We consider asymptotics as the numbers of vertices tend to infinity, and 
thus we assume throughout the paper that we are given, for each n, a se- 
quence d^*^-* = (dj- )i with X^j dj- even. (As usual, we could somewhat 
more generally assume that we are given a sequence n^ ^^ oo and for each 
v a sequence d''^-' = {d^ Yi" ■) For notational simplicity we will usually not 
show the dependency on n explicitly; we thus write d and di, and simi- 
larly for other (deterministic or random) quantities introduced below. All 
unspecified limits and other asymptotic statements are for n ^ cxd. For 
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2 SVANTE JANSON 

example, w.h.p. (with high probabihty) means 'with probabihty tending to 
1 as n — > oo', and — > means 'convergence in probabihty as n — > oo'. Sim- 
ilarly, we use Op and Op in the standard way, always implying n — > oo. For 
example, if X is a parameter of the random graph, X = Op{n) means that 
F{X > en) — > as n — > oo for every e > 0; equivalently, X/n — > 0. 

We may obtain G{n, d) by conditioning the multigraph G*{n, d) on being 
a (simple) graph, i.e., on not having any multiple edges or loops. By Janson 
[9( (with earlier partial results by many authors), 

liminf P(G*(n,d) is simple) > ^^ X] ^? = '^ ( IZ^M ' (^•^) 

i=l \i=l / 

In this case, many results transfer immediately from G*{n, d) to G{n, d), for 
example, every result of the type P(<fn) — *■ for some events £n, and thus 
every result saying that some parameter converges in probability to some 
non-random value. This includes every result in the present paper. 

We will in this paper study the random multigraph G*{n,d); the reader 
can think of doing this either for its own sake or as a tool for studying 
G(n, d). We leave the statement of corollaries for G(n,d), using (jl.ip . to 
the reader. Moreover, the results for G{n, d) extend to some other random 
graph models too, in particular G{n,p) with p ~ X/n and G{n,ni) with 
ni ~ An/2 with A > 0, by the standard device of conditioning on the degree 
sequence; again we omit the details and refer to [10|; [ill; II4I where this 
method is used. 

We will consider percolation of these random (multi)graphs, where we 
first generate a random graph G*{n,d) and then delete either vertices or 
edges at random. (Prom now on, we simply write 'graph' for 'multigraph'.) 
The methods below can be combined to treat the case of random deletion of 
both vertices and edges, which is studied by other methods in e.g. Britton, 
Janson and Martin-Lof [4], but we leave this to the reader. 

To be precise, we consider the following two constructions, given any 
graph G and a probability vr G [0, 1] . 

Site percolation: Randomly delete each vertex (together with all in- 
cident edges) with probability 1 — vr, independently of all other ver- 
tices. We denote the resulting random graph by Gtt.v 
Bond percolation: Randomly delete each edge with probability l—vr, 
independently of all other edges. (All vertices are left.) We denote 
the resulting random graph by Gyr.e- 
Thus TT denotes the probability to be kept in the percolation model. When, 
as in our case, the original graph G itself is random, it is further assumed 
that we first sample G and then proceed as above, conditionally on G. 

The cases vr = 0, 1 are trivial: Gi^ = Gi_e = G, while Go,v = Oj the totally 
empty graph with no vertices and no edges, and Go,e is the empty graph 
with the same edge set as G but no edges. We will thus mainly consider 
0< vr< 1. 
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We may generalize the site percolation model by letting the probability 
depend on the degree of the vertex. Thus, if tt = {TTd)o^ is a given sequence 
of probabilities tt^ E [0, 1], let Gt^.v be the random graph obtained by delet- 
ing vertices independently of each other, with vertex v ^ G deleted with 
probability 1 — 7r(^(^) where d{v) is the degree of v in G. 

For simplicity and in order to concentrate on the main ideas, we will in 
this paper consider only the case when the probability vr (or the sequence 
tt) is fixed and thus does not depend on n, with the exception of a few 
remarks where we briefly indicate how the method can be used also for a 
more detailed study of thresholds. 

The present paper is inspired by Fountoulakis [7[], and we follow his idea of 
deriving results for the percolation models G*(n,d),r,v and G*(n,d)^^e from 
results for the model G*(n,d) without deletions, but for different degree 
sequences d. We will, however, use another method to do this, which we 
find simpler. 

Fountoulakis [7] shows that for both site and bond percolation on G*{n, d), 
if we condition the resulting random graph on its degree sequence d', and 
let n' be the number of its vertices, then the graph has the distribution of 
G*(n',d'), the random graph with this degree sequence constructed by the 
configuration model. He then proceeds to calculate the distributions of the 
degree sequence d' for the two percolation models and finally applies known 
results to G*{n',d'). 

Our method is a version of this, where we do the deletions in two steps. 
For site percolation, instead of deleting a vertex, let us first explode it by 
replacing it by d new vertices of degree 1, where d is its degree; we further 
colour the new vertices red. Then clean up by removing all red vertices. Note 
that the (random) explosions change the number of vertices, but not the 
number of half-edges. Moreover, given the set of explosions, there is a one-to- 
one correspondence between configurations before and after the explosions, 
and thus, if we condition on the new degree sequence, the exploded graph 
is still described by the configuration model. Furthermore, by symmetry, 
when removing the red vertices, all vertices of degree 1 are equivalent, so we 
may just as well remove the right number of vertices of degree 1, but choose 
them uniformly at random. Hence, we can obtain G*{n,d).,^^^ as follows: 

Site percolation: For each vertex i, replace it with probability 1 — vr 
by di new vertices of degree 1 (independently of all other vertices). 
Let d^r^v be the resulting (random) degree sequence, let n be its 
length (the number of vertices), and let n+ be the number of new 
vertices. Construct the random graph G*{h,d.,^^y). Finish by delet- 
ing n+ randomly chosen vertices of degree 1 . 

The more general case when we are given a sequence tt = (Tr^;)^ is handled 
in the same way: 

Site percolation, general: For each vertex i, replace it with proba- 
bility 1 — VTrf. by di new vertices of degree 1. Let d,r,v be the resulting 
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(random) degree sequence, let n be its length (the number of ver- 
tices), and let n+ be the number of new vertices. Construct the 
random graph G*(n,d7r,v)- Finish by deleting n+ randomly chosen 
vertices of degree 1. 

Remark 1.1. We have here assumed that vertices are deleted at random, 
independently of each other. This is not essential for our method, which may 
be further extended to the case when we remove a set of vertices determined 
by any random procedure that is independent of the edges in G{n, d) (but 
may depend on the vertex degrees). For example, we may remove a fixed 
number m of vertices, chosen uniformly at random. It is easily seen that if 
m/n -^ TT, the results of Subsection 12.11 below still hold (with all ttj = vr), 
and thus the results of the later sections hold too. Another, deterministic, 
example is to remove the first m vertices. 

For bond percolation, we instead explode each half-edge with probability 
1 — ^/^^, independently of all other half-edges; to explode a half-edge means 
that we disconnect it from its vertex and transfer it to a new, red vertex of 
degree 1. Again this does not change the number of half-edges, and there 
is a one-to-one correspondence between configurations before and after the 
explosions. We finish by removing all red vertices and their incident edges. 
Since an edge consists of two half-edges, and each survives with probability 
^/^T, this gives the bond percolation model G*(n,d)7r,e where edges are kept 
with probability vr. This yields the following recipe: 

Bond percolation: Replace the degrees di in the sequence d by inde- 
pendent random degrees di ~ Bi((ij, \/n). Add n+ := X^^=i(di — di) 
new degrees 1 to the sequence {di)^, and let d-j^^e be the resulting 
degree sequence and n = n + n+ its length. Construct the random 
graph G*(n, d^r^e)- Finish by deleting n+ randomly chosen vertices 
of degree 1. 

In both cases, we have reduced the problem to a simple (random) modifi- 
cation of the degree sequence, plus a random removal of a set of vertices of 
degree 1. The latter is often more or less trivial to handle, see the applica- 
tions below. We continue to call the removed vertices red when convenient. 

Of course, to use this method, it is essential to find the degree sequence d 
after the explosions. We study this in Section El We then apply this method 
to three different problems: 

Existence of a giant component in the percolated graph, i.e., what is called 
percolation in random graph theory (Section (S]). Our results include and 
extend earlier work by Fountoulakis [TJ, which inspired the present study, 
and some of the results by Britton, Janson and Martin-Lof [J]. 

Existence of a k-core in the percolated graph (Section Hj). We obtain a 
general result analogous to (and extending) the well-known result by Pittel, 
Spencer and Wormald p^] for G{n,p). We study the phase transitions that 
may occur in some detail and show by examples that it is possible to have 
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several, and even an infinite number of, different phase transitions as the 
probability vr increases from to 1. 

Bootstrap percolation in random regular graphs (Section [5]), where we 
obtain a new and simpler proof of results by Balogh and Pittel [l| . 

For a graph G, let v{G) and e(G) denote the numbers of vertices and 
edges in G, respectively, and let Vj{G) be the number of vertices of degree j, 
j > 0. We sometimes use G*{n, d),^ to denote any of the percolation models 
G*(n,d)^,v, G*(n,d)„,v or G*(n,d)^,e- 

2. The degree sequence after explosions 

Let rij := #{i < n : di = j}, the number Vj{G*{n,d)) of vertices of degree 
j in G*(n,d). Thus Y27Lo''^j ~ "-• ^^ assume for simplicity the following 
regularity condition. 

Condition 2.1. There exists a probability distribution {pj)JLq with finite 
positive mean A := YljJPj ^ (^1) oo) such that (as n -^ oo) 

Uj/n^Pj, j > 0, (2.1) 

and 

Eoo 
7=0 J"i 



n 



Y.^Pr (2.2) 

3=0 



Note that, in order to avoid trivialities, we assume that A > 0, which is 
equivalent to po < 1- Thus, there is a positive fraction of vertices of degree 
at least 1. 

Note that "^jjrij = Yli^i equals twice the number of edges in G*(n, d), 
and that ()2.2p says that the average degree in G*{n,d) converges to A. 

Let the random variable D = Dn be the degree of a random vertex in 
G*(re, d), thus Z)„, has the distribution P(-D„ = j) = rij/n, and let D be 
a random variable with the distribution (pj)(f . Then (j2.ip is equivalent 

to Dn -^ D, and ([22]) is Ei)„ ^ A = ED. Further, assuming (f2T|l . 



(j2.2p is equivalent to uniform integr ability of Dn, or equivalently uniform 
summability (as n — > cxd) of "^j jnj/n, see for example Gut [g. Theorem 
5.5.9 and Remark 5.5.4]. 

Remark 2.2. The uniform summability of ^ • jrij/n is easily seen to imply 
that if ^ is any (random or deterministic) subgraph on G*{n, d) with v{H) = 
o{n), then e{H) = o{n), and similarly with Op{n). 

We will also use the probability generating function of the asymptotic 
degree distribution D: 

oo 

gD{x):=Ex^ = Y,Pjx', (2-3) 

i=o 

defined at least for Ixl < 1. 
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We perform either site or bond percolation as in Section [H by the explo- 
sion method described there, and let hj := ^{i < n : di = j} be the number 
of vertices of degree j after the explosions. Thus 

oo 

^nj=n. (2.4) 

i=o 

It is easy to find the distribution of (hj) and its asymptotics for our two 
percolation models. 

2.1. Site percolation. We treat the general version with a sequence tt. 
Let n^ be the number of vertices of degree j that are not exploded. Then 

n° ~ Bi(nj,7rj) (independent of each other), (2-5) 

oo 

n+ = ^j{nj-n°), (2.6) 

n,=n°, j^l, (2.7) 

hi = nl + n+. (2.8) 

By the law of large numbers, n", = tijitj + Op{n) and thus, using the assump- 
tion ()2.ip and the uniform summability of ^jjnj/n (which enables us to 
treat the infinite sums in (I2.10p and ()2.13p by a standard argument), 

n° = rijTTj + Op{n) = ttjPjTI + Op{n), (2.9) 

oo 

n+ = ^i(l - T^j)Pjn + Op{n), (2.10) 

j=0 
"j = ^jPj^ + Op{n), j 7^ 1, (2.11) 

oo 

hi = (nipi + ^ j(l - iTj)pjjn + Op{n), (2.12) 

j=0 

oo 

n = X](^J +-^(^ ~ '^j))Pj^ + Opin). (2.13) 

j=0 



We can write (|2.13p as 



n p 
n 






Further, by (pTTT) and (|2T2]) . 

n '^'- lc-^(vriPi+Er=iJ'(l-vr,)p,), j = I. ^ ^ 
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Since hj < rij for j > 2 and h > n — uq, the uniform suniniability of juj/n 
implies uniform summability of jhj/h, and thus also 

2^7=0 J "-j p 



n 

3=0 



X:=Y.3Pj<°^- (2-16) 



Hence Condition 12.11 holds . in probability, for the random degree sequence d 
too. Further, the total number of half-edges is not changed by the explosions, 
and thus also, by ([2311) and dZ^J, 



j=oJnj _ Ej=oJrij _ n Ej=oJ"j p ^ ^-i x_. (2.17) 

h n h n 

hence (or by (|2.15|) ). 

A = C^X. (2.18) 

In the proofs below it will be convenient to assume that (j2.15p and (j2.16p 
hold a.s., and not just in probability, so that Condition 12. II a. s. holds for d; 
we can assume this without loss of generality by the Skorohod coupling the- 
orem [13 . Theorem 4.30]. (Alternatively, one can argue by selecting suitable 
subsequences.) 

Let D have the probability distribution {pj), and let qj^ be its probability 
generating function. Then, by (j2.15p . 

oo oo oo oo 

j=0 j=0 j=0 j=0 

(2.19) 
In particular, if all ttj = n, 

Cgf){x)=TrgDix) + {l-7r)Xx, (2.20) 

where now (^ = tt + (1 — tt)X. 

2.2. Bond percolation. For bond percolation, we have explosions that do 
not destroy the vertices, but they may reduce their degrees. Let hji be the 
number of vertices that had degree I before the explosions and j after. Thus 
fij = Yli>j ^ji fo'^ J 7^ 1 s-iid hi = Ylii>i ^11 -|- n+. A vertex of degree I will 
after the explosions have a degree with the binomial distribution Bi(l, tt^'^), 
and thus the probability that it will become a vertex of degree j is the 
binomial probability bij{'iT^''^), where we define 

bijip) := P(Bi(/,p) = j) = Qpi{l-py-^. (2.21) 

Since explosions at different vertices occur independently, this means that, 
for l>j>0, 

hji ~ Bi(n;,6;j(7r^/^)) 
and thus, by the law of large numbers and (j2.1|) . 

hji = bij{TT^^'^)pin + Op{n). 
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Further, the number n^- of new vertices equals the number of explosions, 
and thus has the binomial distribution Bi(^j/ni, 1 — vr^'^). Consequently, 
using also (|2.2p and the uniform summability of X^j i^i/''^> 

n+ = Y^ lni{l - TT^/^) + o^(„) = (1 - 7r^^^)Xn + Op{n), (2.22) 

I 

^j = Yl "j' = Yl bij{Tr^^'')pin + Op{n), j ^ 1, (2.23) 

l>j l>j 

ni = ^ nii + n+ = ^ bn{TT^/^)pin + (l - TT^/^)Xn + Op{n), (2.24) 

n = n + n+ = n+(l- 7r^/^)An + Op(n). (2.25) 

In analogy with site percolation we thus have, by (j2.25p . 
h 
n 
and further, by (12:231) and OTMl) . 



-^(:=l + (l-vri/2)A (2.26) 



n '^^■"\C-^(Em^n(^^/>/ + (l--^/^)A), j = l. ^'-"'^ 

Again, the uniform summability of juj/n implies uniform summability of 
jfij/h, and the total number of half-edges is not changed; thus ()2.16p . ()2.17p 
and (j2.18p hold, now with (^ given by (|2.26p . Hence Condition 12.11 holds in 
probability for the degree sequences d in bond percolation too, and by the 
Skorohod coupling theorem we may assume that it holds a.s. 

The formula for pj is a bit complicated, but there is a simple formula for 
the probability generating function gj=). We have, by the binomial theorem, 
Ej</ hj{T^)x^ = (1 - vr + vrx)', and thus ^^TT!) yields 

oo 

Cg^{x) = ^(1 - ^1/2 + ^i/2x)V/ + (1 - vri/2)Ax 

1=0 

= goil - ^^/^ + ^^'"^x) + (1 - ^^''^)\x. (2.28) 

3. Giant component 
The question of existence of a giant component in G(n,d) and G*(n, d) 



was answered by Molloy and Reed 15[], who showed that (under some weak 



technical assumptions) a giant component exists w.h.p. if and only if (in the 
notation above) ¥.D{D — 2) > 0. (The term giant component is in this paper 
used, somewhat informally, for a component containing at least a fractions 
e of all vertices, for some small e > that does not depend on n.) They 
further gave a formula for the size of this giant component in Molloy and 
Reed |16l] . We will use the following version of their result, given by Janson 



and Luczak 1^, Theorem 2.3 and Remark 2.6]. Let, for any graph G, Ck{G) 
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denote the k:th largest component of G. (Break ties by any rule. If there 
are fewer that k components, let C^ := 0.) 

Proposition 3.1 ([16|;ll4l)' Consider G*{n,d), assuming that Condition \2.1\ 
holds and pi > 0. Let C^ ■= Ck{G*{n,d)) and let goix) be the probability 
generating function in (j2.3|) . 

(i) // E D{D - 2) = Y.J JiJ - '^)Pj > 0, then there is a unique ( G (0, 1) 
such that 5'^(^) = A^, and 

viCi)/n^l-gD{0>0, (3.1) 

Vj{Ci)/n -^ pj{l - ^^), for every j > 0, (3.2) 

e(Ci)/n^iA(l-e'). (3.3) 

Furthermore, v{C2)/n — > and e{C'i)ln — > 0. 
(ii) // ED(Z) - 2) = XljiO' - 2)Pj < 0, then v(Ci)/n -^ and 
e{Ci)/n ^ 0. 

Remark 3.2. ED^ = cxo is allowed in Proposition I3.]] r i)[ 



Remark 3.3. In Proposition 13. l|(n)| where 'ED{D — 2) < and pi < 0, for 
0< x< 1 



Ax - g'uix) = ^jpj{x -x^ ^) = pi{x - 1) +x^jpj{l - x 



3=1 i=2 



<Pi{x-l)+xY^jPj{j-2){l 



i=2 

oo 

< ^3{i - 2)pjx{l - x) =ED{D - 2)x(l - x) < 0. 
i=i 

Hence, in this case the only solution in [0, 1] to g'^iC) = A^ is ^ = 1, which 
we may take as the definition in this case. 

Remark 3.4. Let D* be a random variable with the distribution 

nD* = j) = (j + 1) nD = j + 1)/A, j > 0; 

this is the size-biased distribution of D shifted by 1, and it has a well-known 
natural interpretation as follows. Pick a random half-edge; then the number 
of remaining half-edges at its endpoint has asymptotically the distribution of 
D*. Therefore, the natural (Galton- Watson) branching process approxima- 
tion of the exploration of the successive neighbourhoods of a given vertex is 
the branching process X with offspring distributed as D*, but starting with 
an initial distribution given by D. Since 
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the equation g'^i^) = A^ in Proposition l3.]|Fi)] can be written gD*{0 = Ci 
which shows that (, has an interpretation as the extinction probabihty of 
the branching process X with offspring distribution D*, now starting with a 
single individuah (This also agrees with the definition in Remark 13.31 for the 
case Proposition I3.1|ni)| ) Thus goiO ^^ (|3.ip is the extinction probability 
of X. Note also that 

^j^.^^D{D-l) _EDiD-l) 



A ED 

so the condition ED{D - 2) > 0, or equivalently ED{D - 1) > EL>, is 
equivalent to ED* > 1, the classical condition for the branching process to 
be supercritical and thus have a positive survival probability. 

The intuition behind the branching process approximation of the local 
structure of a random graph at a given vertex is that an infinite approximat- 
ing branching process corresponds to the vertex being in a giant component. 
This intuition agrees also with the formulas (|3.2p and (j3.3p . which refiect 
the fact that a vertex of degree j [an edge] belongs to the giant component 
if and only if one of its j attached half-edges [one of its two constituent 
half-edges] connects to the giant component. (It is rather easy to base rig- 
orous proofs on the branching process approximation, see e.g. [J], but in the 
present paper we will only use the branching process heuristically.) 

Consider one of our percolation models G*{n, d),^, and construct it using 
explosions and an intermediate random graph G*{h,d) as described in the 
introduction. (Recall that d is random, while d and the limiting probabilities 
Pj and pj are not.) Let Cj := Cj (G*(n, d),^) and Cj := Cj(G*(n, d)) denote 
the components of G*(n, d)^, and G*(n, d), respectively. 

As remarked in Section [2l we may assume that G*(n,d) too satisfies 
Condition l2.1l with pj replaced by pj. (At least a.s.; recall that d is random.) 
Hence, assuming pi > 0, if we first condition on d, then Proposition 13.11 
applies immediately to the exploded graph G*{n, d). We also have to remove 
n+ randomly chosen "red" vertices of degree 1, but luckily this will not 
break up any component. Consequently, if ED(D — 2) > 0, then G*(fi, d) 
w.h.p. has a giant component Ci, with v{Ci), Vj{Ci) and e(Ci) given by 
Proposition [3]T] (with pj replaced by pj), and after removing the red vertices, 
the remainder of Ci is still connected and forms a component C in G*{n, d)-,^. 
Furthermore, since ED{D — 2) > 0, pj > for at least one j > 2, and it 
follows by (j3.2p that Ci contains en + Op{n) vertices of degree j, for some 
c > 0; all these belong to C (although possibly with smaller degrees), so C 
contains w.h.p. at least cn/2 vertices. Moreover, all other components of 
G*(n, d)^ are contained in components of G*(n,d) different from Ci, and 
thus at most as large as C2, which by Proposition 13.11 has Op{n) = Op{n) 
vertices. Hence, w.h.p. C is the largest component Ci of G*(n, d)^, and this 
is the unique giant component in G*{n,d)n- 
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Since we remove a fraction rij^/hi of all vertices of degree 1, we remove 
by the law of large numbers (for a hypergeometric distribution) about the 
same fraction of the vertices of degree 1 in the giant component Ci. More 
precisely, by ()3.2I1 . Ci contains about a fraction 1 — ^ of all vertices of degree 
1, where g'f.iO = -^Ci hence the number of red vertices removed from Ci is 

(l-C)n+ + Op(n). (3.4) 

By dSlD and SSlh . 

v{Ci) = v{Ci)-{l-^)n+ + Op{n) = n(l-5(^(^)) -n+ + n+^ + Op{n). (3.5) 

Similarly, by (j3.3p and (j3.4p . since each red vertex that is removed from Ci 
also removes one edge with it, 

e(Ci) = e(Ci) - (1 - e)n+ + Op{n) = fXh{l - ^^) - (1 - On+ + Op{n). (3.6) 

The case E,D{D — 2) < is even simpler; since the largest component 
Ci is contained in some component Cj of G*{h,d), it follows that v{Ci) < 
v{Cj) < v{Ci) = Op{h) = Op{n). 

This leads to the following results, where we treat site and bond percola- 
tion separately and add formulas for the asymptotic size of Ci . 

Theorem 3.5. Consider the site percolation model G*{n, d)7r,v) <ind suppose 
that Condition \2.1\ holds and that n = (7rrf)[^ with < vr^ < 1; suppose 
further that there exists j > 1 such that pj > and ttj < 1. Then there is 
w.h.p. a giant component if and only if 

oo oo 

i=0 j=0 

(i) // (|3.7p holds, then there is a unique ^ = Cvi"^) £ (0, 1) such that 

OD 

Y^jn,p,{l-e-^) = X{l-0 (3.8) 

i=i 

and then 

oo 

v{Ci)/n ^ xv{7v) := ^^,p,(l - e^) > 0, (3.9) 

e{C,)/n -^ ^v(7r) := (1 - E^'^^'^^' " -^ Y.^Pr (3-10) 



2 

Furthermore, v{C2)/n — > and e[C'i)ln — > 0. 
(ii) // ()3.7p does not hold, then v[Ci)/n — > and e{Ci)/n — > 0. 
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Proof. We apply Proposition 13.11 to G*{h,d) as discussed above. Note that 
pi > by (I2.15P and the assumption (1 — itj)pj > for some j. By (j2.15p . 

oo oo oo 

CEDiD-2) = C^jij - 2)pj = Y,J{3 - '^)^jPj - E-?'(l - ^i)Pi 
i=o j=i j=i 

oo oo 

Hence, the condition E,D{D — 2) > is equivalent to ([3 



In particular, it follows that f (C2) = Op{n) in (i) and v{Ci) = Op{n) in 



(ii) That also e(C2) = Op{n) in (i) and e(Ci) = Op{n) in (ii) follows by 



Remark [2^2] applied to G*(n, d). 

It remains only to verify the formulas (|3.8p - (|3.10p . The equation g'f,{0 = 
A^ is by (|2.18p equivalent to Cg'f^{^) = A^, which can be written as (|3.8p by 
()2.15p and a simple calculation. 

By dSSI), using (IXTU|) . (IXTip and (1X1^1) . 

00 

v{Ci)/n ^ c - CgdO - (1 - E-?'(i - ^j)^i 

00 00 CO 

i=o i=o j=0 

00 

Similarly, by ([32]), (|2T8D . (|2TiD and (|230]) . 

00 
e(Ci)/n ^ iA(l - e') - (1 - 5^ j(l - vr,)p, 

i=i 

In the standard case when all Hd = vr, this leads to a simple criterion, 
which earlier has been shown by Britton, Janson and Martin-Lof [4] and 
Fountoulakis [7] by different methods. (A modification of the usual branch- 
ing process argument for G*(n, d) in [4] and a method similar to ours in 

00 

Corollary 3.6 ([1; 0]). Suppose that Condition \2J\ holds and < vr < 1. 
Then there exists w.h.p. a giant component in G*(n,d)7r,v if o,''^d only if 

ED 

TT > TTc := z r. D (3.11) 

ED{D-1) ^ ' 
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Remark 3.7. Note that tTc = is possible; this happens if and only if 
KD'^ = cxD. (Recall that we assume < ED < cxd, see Condition 12.11 ) Fur- 
ther, TTc > 1 is possible too: in this case there is w.h.p. no giant component 
in G*{n, d) (except possibly in the special case when pj = for all j ^ 0, 2), 
and consequently none in the subgraph G*{n,d)-,^. 

Note that by ^JB), tTc € (0, 1) if and only if ED < ED{D - 1) < oo, 
i.e., if and only if < ED{D - 2) < oo. 



Remark 3.8. Another case treated in [J] (there called El) is tt^ = a'^ for 
some a G (0, 1). Theorem 13.51 gives a new proof that then there is a giant 
component if and only if X^^x j(j — l)a^Pj > A, which also can be written 



a 



2„n 



g'l){a) > X = g'jj{l). (The cases E2 and A in [J] are more complicated and 
do not follow from the results in the present paper.) 

For edge percolation we similarly have the following; this too has been 
shown by Britton, Janson and Martin-Lof [4] and Fountoulakis [7\. Note 
that the permutation threshold vr is the same for site and bond percolation, 
as observed by Fountoulakis [Tj. 

Theorem 3.9 ([J;l7|). Consider the bond percolation model G*{n,d)T^^e, and 
suppose that Condition \2. 1\ holds and that < vr < 1. Then there is w.h.p. 
a giant component if and only if 

ED 

TT > TTc := z r. (3.12) 

ED(D-l) ^ ' 

(i) If (|3.12p holds, then there is a unique S, = (.ei^^) G (0, 1) such that 
n'/'g'j, (1 - 7rV2 + ^y^^) + (i _ vr^^) a = A^, (3.13) 

and then 

v{Ci)/n ^ xe(vr) := 1 - g^ (l - ^'/' + ^'^'e) > 0, (3.14) 

e(Ci)/n ^ ^e(vr) := tt'^Hi - e)A - iA(l - O'- (3-15) 

Furthermore, v{C2)/n — > and e{C2)/n — > 0. 
(ii) If (|3.12p does not hold, then v{Ci)/n — > and e{Ci)/n — > 0. 

Proof. We argue as in the proof of Theorem 13.51 noting that pi > by 
(12:2711 . By (ESSD, 

CED(D - 2) = Cglii) - Cg'f,{i) = TrgUi) - ^'^'g'Dii) - (i - ^'/')A 

= 7rED(D-l)- A, 

which yields the criterion (13.120 . Further, if ()3.12p holds, then the equation 
g'f){C) = A.^, which by ()2.18p is equivalent to C5'rs(0 = CA^ = A.^, becomes 

^1^ by (f2:28]l. 

By ^, (|2:26D . (12:221) and (12:281) . 

v{C,)/n ^ C - CgoiO - (1 - 0(1 - ^'/')A = 1 - 5D (1 - vrV2 + vr^s^)) , 
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which is (IXTill . Similarly, I^Ml, dZ^S]), (|2J8]l and (12:221) yield 

e(Ci)/n ^ iA(l - e) - (1 - e)(l - vri/')A = vri/2(l - ^A - iA(l - O', 
which is (j3.15p . The rest is as above. D 

Remark 3.10. It may come as a surprise that we have the same criterion 
(j3.1ip and (j3.12p for site and bond percolation, since the proofs above arrive 
at this equation in somewhat different ways. However, remember that all re- 
sults here are consistent with the standard branching process approximation 
in Remark 13.41 (even if our proofs use different arguments) and it is obvious 
that both site and bond percolation affect the mean number of offspring in 
the branching process in the same way, namely by multiplication by vr. Cf. 
[J], where the proofs are based on such branching process approximations. 

Define 

Pv = Pv(7r) := 1 - Cv(7r) and pe = Pe(vr) := 1 - Ce(vr); (3.16) 

recall from Remark 13.41 that ,^v and .^e are the extinction probabilities in the 
two branching processes defined by the site and bond percolation models, 
and thus py and pe are the corresponding survival probabilities. For bond 
percolation, (I3.13p - (l3.15p can be written in the somewhat simpler forms 

vr^/'5i,(l-vrVV)=A(vrV2_p,), (3.17) 

v{Ci)/n ^ Xe{TT) := 1 - 5D (1 - vri/Ve(^)) , (3.18) 

e(Ci)/n ^ pe(^) := 7^^^^Xpe{7T) - IXpeinf. (3.19) 

Note further that if we consider site percolation with all ttj = vr, (j3.8p can 
be written 

7r(A-ff^(l-/>v)) =Apv (3.20) 

and it follows by comparison with (j3.17p that 

p,{n) = 7ri/Ve(vr). (3.21) 

Furthermore, ([SJD, (|310]) . (|3T8D and (|3T9D now yield 

Xv(vr) = 7r(l - goiU^))) = 7t{1 - <7d(1 - Pv(vr))) = ^Xe(^), (3.22) 
Pv(vr) = 7rApv(vr) - ^Apv(vr)^ = vrpe(7r). (3.23) 

We next consider how the various parameters above depend on tt, for both 
site percolation and bond percolation, where for site percolation we in the 
remainder of this section consider only the case when all ttj = vr. 

We have so far defined the parameters for vr G (tTc, 1) only; we extend 
the definitions by letting ^v '■= Ce '■= 1 and p^ '■= Pe '■= Xv '■= Xe '■= pv '■= 
Pe '■= for vr < TTc, noting that this is compatible with the branching process 
interpretation of ^ and p in Remark l3.4l and that the equalities in (J3.8p ~ (j3.23p 
hold trivially. 
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Theorem 3.11. Assume Condition [Kj^ The functions ^y, py,Xv, fJ-v^Ce, Pe, 
Xe, pe are continuous functions of it £ (0, 1) and are analytic except at it = 
TTc. (Hence, the functions are analytic in (0,1) if and only if tTc = or 

TTc > I.) 

Proof. It suffices to sliow tliis for ^v! the result for tlie otlier functions then 
follows by (j3.16p and (J3.2ip - (j3.23p . Since the case tt < tTc is trivial, it suffices 
to consider vr > tTc, and we may thus assume that < vTc < 1. 

If vr e (tTc,!), then, as shown above, g'f,{(,\/) = ^S,v, or, equivalently, 
Gi(,v,Tr) = 0, where G(^,7r) := g'fyi(,)/(,~^ is an analytic function of (^,vr) G 
(0, 1)^. Moreover, G(^,7r) is a strictly convex function of ^ G (0, 1] for any 
vr G (0, 1), and G(^v, vr) = G(l, vr) = 0; hence 9G(^v, vr)/a^ < 0. The implicit 
function theorem now shows that ^vC^r) is analytic for vr G (tTc, 1). 

For continuity at tTc, suppose tTc G (0,1) and let ^ = liuin~*ooS,v{'^n) 
for some sequence 7r„ — > tt^- Then, writing D{7r) and A(7r) to show the 
dependence on n, g'- -.Hvi'^n)) = ^{'^n)Cv{'^n) and thus by continuity, e.g. 

using (12:28]) . g'-. .{i) = A(7rc)|. However, for vr < vTc, we have Ei)(Z)-2) < 

and then ^ = 1 is the only solution in (0, 1] of g'f.{0 — -^Ci hence ^ = 1. 
This shows that ^v(7r) ^ 1 as vr ^ vTc, i.e., ^v is continuous at vTc. □ 

Remark 3.12. Alternatively, the continuity of ^v in (0, 1) follows by Re- 
mark 13.41 and continuity of the extinction probability as the offspring dis- 
tribution varies, cf. [J, Lemma 4.1]. Furthermore, by the same arguments, 
the parameters are continuous also at vr = and, except in the case when 
Po + P2 = 1 (and thus D = 1 a.s.), at vr = 1 too. 

At the threshold vTc, we have linear growth of the size of the giant com- 
ponent for (slightly) larger vr, provided KD"^ < oo, and thus a jump discon- 
tinuity in the derivative of ^vjXvi • • • • More precisely, the following holds. 
We are here only interested in the case < vTc < 1, which is equivalent to 
< KD{D - 2) < oo, see Remark [321 

Theorem 3.13. Suppose that < ED(D - 2) < oo; thus < tt^ < 1. If 
further KD"^ < oo, then as e \0, 

2ED{D-1) 2{ED{D-l))'^ 

Pv(vrc + e) ~ ^^jE^p _ ^)(^ -2f-ED -EDiD - 1){D - 2)^ ^^ ' 

> . ^ 2ED-ED(D-l) 

Xv(vrc + e) ~ /iv(vrc + e) ~ 7rcApv(^c + e) ~ eD(D - l)(D - 2f ' ^ ' 

Similar results for p^, Xe, Me follow by ()3.2ip - ()3.23p . 

Proof For vr = vrc -F e \ vrc, by g^(l) =ED{D - 1) = A/vrc, see (l3TT]l . and 
dS^QD, 

e5'D(l)Pv = (vr-vrc)5f^(l)pv = vrs-^(l)/3v-A/3v = vr(5-^(l)pv-A-F5^(l-pv))- 

(3.26) 
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Since KD^ < oo, qd is three times continuously differentiable on [0, 1], and 
a Taylo r exp ansion yields g^(l - p^) = X - pyg'^{l) + p^5^'(l)/2 + o{pl). 
Hence, (I3.26P yields, since pv > 0, 

egW^) = 7rpv9D(l)/2 + o{p,) = 7rcPv<7D(l)/2 + o(pv). 

Thus, noting that 5^(1) =¥.D{D - 1) and c/^'(l) = ¥.D{D - 1){D - 2) > 
(sinceEL>(D-2) > 0), 

2g^(l)_ 2ED{D-1) 



^^ TT.g'l^il) 7r,EDiD-l){D-2) 
which yields dS^D- Finally, ([3:25]) follows easily by ([3:22]) and ([3:23]) . D 

If KD^ = oo, we find in the same way a slower growth of Pvi"^), X'ji'^)i 
Pv{tt) at TTc- As an example, we consider D with a power law tail, pk ~ ck~'^, 
where we take 3 < 7 < 4 so that ED^ < 00 but ED^ = 00. 

Theorem 3.14. Suppose that pk ~ ck~'^ as k ^ 00, where 3 < 7 < 4 and 
c > 0. Assume further that KD(D — 2) > 0. Then tt^ £ (0, 1) and, as e \0, 

Ei?(Z)-l)V/^"-'^i/(,-3) 



^^^"^+"^~^vr.r(2-7)' 

Xv(vrc + e) ~ Aiv(vrc + e) ~ ttcApvIttc + e) 

Similar results for p^, Xe; /"e follow by ([3.2ip - ([3.23p . 

Proof. We have, for example by comparison with the Taylor expansion of 

(i-(i-t)r', 

00 

ffD(l - t) = 5] A:(A: - 1)(A; - 2)^^(1 - t)'= ~ cr(4 - j)t^-\ t \ 0, 

A;=3 

and thus by integration 

5^(1) - g'Ul - t) ~ cr(4 - 7)(7 - 3)-4^-3 = c|r(3 - 7)|t^-^ 
and, integrating once more, 

pWM -{X-g'D^- Pv)) ~ cr(2 - 7)/or'- 

Hence, ([3.26P yields 

e5D(lK~c7rcr(2-7)pr2 

and the results follow, again using ([3.22p and (|3.23p . D 
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4. /c-CORE 

Let A; > 2 be a fixed integer. The fc-core of a grapli G, denoted by 
Corefc(G), is tlie largest induced subgraph of G with minimum vertex de- 
gree at least k. (Note that the /c-core may be empty.) The question whether 
a non-empty fe-core exists in a random graph has attracted a lot of atten- 
tion for various models of random graphs since the pionee ring papers by 
Bollobas [2], Luczak [iJ] and Pittel, Spencer and Wormald [171] for G{n,p) 
and G{n,m); in particular, the case of G(n, d) and G*{n,d) with given de- 
gree sequences have been studied by several authors, see Janson and Luczak 
[id . [ll[ and the references given there. 

We study the percolated G*{n,d)T^ by the exposion method presented in 
Section [TJ For the A;-core, the cleaning up stage is trivial: by definition, 
the A:-core of G*(n, d) does not contain any vertices of degree 1, so it is 
unaffected by the removal of all red vertices, and thus 

Corefc(G*(n,d)^) = Corefc(G*(n,d)). (4.1) 

Let, for < p < 1, Dp be the thinning of D obtained by taking D points 
and then randomly and independently keeping each of them with probability 
p. Thus, given D = d, Dp ~ Bi(d, p). Define, recalling the notation (|2.2ip . 

oo oo 

h{p) := E{Dpl[Dp > k]) = Y.Y.jP^h^{p), (4.2) 

j=k l=j 

oo oo 

hi{p) := ¥{Dp > k) = Y^Y.pibij{p). (4.3) 

j=k i=j 

Note that Dp is stochastically increasing in p, and thus both h and hi are 



increasing inp, with h{0) = /ii(0) = 0. Note further that h{l) = YlT=kJPj — 
A and /ii(l) = YlT=kPj — ^^ with strict inequalities unless pj = for all 
j = 1, . . . ,k — 1 ov j = 0,1, . . . ,k — 1, respectively. Moreover, 

fc-i 



h{p) =EDp- E{Dpl[Dp <k-l])=EDp- Y^J^i^p = j) 
= ^P-YYjPi( .JP'(I-P)' 



fc-i 



k-l 



= ^P-T.jj^_9'i\l-P). (4.4) 

Since goiz) is an analytic function in {z : |z| < 1}, (j4.4p shows that h{p) is 
an analytic function in the domain {p : |p — 1| < 1} in the complex plane; 
in particular, h is analytic on (0,1]. (But not necessarily at 0, as seen by 
Example 14.131 ) Similarly, hi is analytic on (0, 1]. 

We will use the following result by Janson and Luczak lOJ, Theorem 2.3]. 
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Proposition 4.1 ([lQ[]). Suppose that Condition \2.1\ holds. Let k > 2 be 
fixed, and let Core^ be the k-core of G*{n,d). Let p := max{p G [0,1] : 
hip)=Xp^}. 

(i) // h{p) < Xp^ for all p G (0, 1], which is equivalent to p = 0, then 
Core^. has Op{n) vertices and Op{n) edges. Furthermore, if also k>3 
and Yl^=i ^"'^^ — 0{n) for some a > 0, then Core^ is empty w.h.p. 
(ii) // h{p) > Xp^ for some p G (0, 1], which is equivalent to p £ (0, 1], 
and further p is not a local maximum point of h{p) — Xp^ , then 

v{CoTel)/n ^ hi{p) > 0, (4.5) 

oo 

vj{Corel)/n -^ F{Dp = j) = Y.pihj{p), j > k, (4.6) 

e{CoTel)/n ^ h{p)/2 = Xf/2. (4.7) 

Remark 4.2. The result (j4.6p is not stated explicitly in [lO|, but as re- 
marked in lll|, Remark 1.8], it follows immediately from the proof in 10[] of 
(14. Sp . (Cf. [5] for the random graph G(n,m).) 



Remark 4.3. The extra condition in (ii) that p is not a local maximum 
point of hip) - Xp' is actually stated somewhat differently in 0, viz. as 
Xp^ < h{p) in some interval {p — £,p). However, since g{p) := h{p) — Xp^ 
is analytic at p, a Taylor expansion at p shows that for some such interval 
(j) — e,p), either g{p) > or g{p) < throughout the interval. Since g{p) = 
and g{p) < for p < p < 1, the two versions of the condition are equivalent. 
The need for this condition is perhaps more clearly seen in the percolation 
setting, cf. Remark |4. 81 



There is a natural interpretation of this result in terms of the branching 
process approximation of the local exploration process, similar to the one 
described for the giant component in Remark 13.41 For the A:-core, this was 
observed already by Pittel, Spencer and Wormald [l7|], but (unlike for the 
giant component), the branching process approximation has so far mainly 
been used heuristically; the technical difficulties to make a rigorous proof 
based on it are formidable, and have so far been overcome only by Riordan 
[13] for a related random graph model. We, as most others, avoid this 
complicated method of proof, and only identify the limits in Proposition 14.11 
(which is proved by other, simpler, methods in |lO|]) with quantities for the 
branching process. Although this idea is not new, we have, for the random 
graphs that we consider, not seen a detailed proof of it in the literature, so 
for completeness we provide one in Appendix lAl 

Remark 4.4. If k = 2, then (lOj) yields 

oo 

Kp) =Xp- ^pilp{l - py-^ = Xp- pg'uil - p) (4.8) 
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and thus 

h{p) -A/ =p[X{l-p) - g'jj(l-p)). 

It follows that p = 1 — ^, where S, is as in Proposition 13.11 and Remark I3.3| 
i.e., by Remark 13.41 p = p, the survival probability of the branching process 
X with offspring distribution D*. (See AppendixRlfor further explanations 
of this.) 

We now easily derive results for the A;-core in the percolation models. For 
simplicity, we consider for site percolation only the case when all vr^ are 
equal; the general case is similar but the explicit formulas are less nice. 

Theorem 4.5. Consider the site percolation model G*(n,d)^^v with < 
IT < 1, and suppose that Condition \2.1\ holds. Let k > 2 be fixed, and let 
Core^ be the k-core o/G*(n, d)7r,v Let 

e^o ■ f V ( Kp) \ (^ a\ 

VTc = vr^ ' := mf -—- = sup —^ . (4.9) 

o<P<i h{p) \o<p<i Ap2 J 

(i) If IT < TTc, then Core^ has Op{n) vertices and Op{n) edges. Further- 
more, if also k > 3 and Y17=i ^°"^^ ~ 0{n) for some a > 0, then 
Core^ is empty w.h.p. 

(ii) //vr > TTc, then w.h.p. Core^ is non-empty. Furthermore, ifp = p{'K) 
is the largest p <1 such that h{p)/{\p^) = vr^^, and p is not a local 
maximum point ofh{p)/{Xp^) in (0,1], then 

v{CoTe%)/n — > ■Khi{p) > 0, 
7;,-(Core^)/n ^ ^P(Dj^ = j), j > k, 

e{CoTel)/n ^ Tih{p)/2 = \f /2. 

Theorem 4.6. Consider the bond percolation model G*{n,A)-,^^^ with < 
vr < 1, and suppose that Condition \2.1\ holds. Let k > 2 be fixed, and let 
Core^ be the k-core of G*{n,d)T^^e. Let tTc = vtc be given by (j4.9p . 

(i) If IT < TTc, then Core^ has Op{n) vertices and Op{n) edges. Further- 
more, if also k > 3 and Y17=i ^ ~ 0{n) for some a > 0, then 
Core^ is empty w.h.p. 
(ii) //vr > TTc, then w.h.p. Core^ is non-empty. Furthermore, ifp = p{'K) 
is the largest p <1 such that h{p)/{\p'^) = vr^^, and p is not a local 
maximum point of h{p) / {Xp^) in (0,1], then 

v^CoieD/n — > hi{p) > 0, 
Vj{Covel)/n^¥{Dp=j), j>k, 
e{CoTel)/n ^ h{p)/2 = Xf/{27T). 

For convenience, we define ip{p) := h{p)/p'^, < p < 1. 
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Remark 4.7. Since h{p) is analytic in (0, 1), there is at most a countable 
number of local maximum points of ip{p) := h{p)/p'^ (except when h{p)/p^ is 
constant), and thus at most a countable number of local maximum values of 
h{j))j]? . Hence, there is at most a countable number of exceptional values of 
TT in part (ii) of Theorems 14.51 and 14. 61 At these exceptional values, we have 
a discontinuity of p(7r) and thus of the relative asymptotic size 7r/ii(p(7r)) 
or hi(p(TT)) of the A;-core; in other words, there is a phase transition of the 
fc-core at each such exceptional n. (See Figures [l] and O) Similarly, if 99 
has an inflection point at ^(vr), i.e., if ip'{p) = ^"{p) = • • • = Lp^^^>{p) = 
and 99'^^+^) (p) < for some £ > 1, then ^(vr) and /ii(p(7r)) are continuous 
but the derivatives of p(7r) and /ii(p(7r)) become infinite at this point, so we 
have a phase transition of a different type. For all other tt > vTc, the implicit 
function theorem shows that p{'k) and /ii(p(7r)) are analytic at vr. 

Say that p is a critical point of ip if (p'{p) = 0, and a had critical point if 
further, p G (0,1), ip{p) > X and ip{p) > ip{p) for all p £ (p, 1). It follows 
that there is a 1-1 correspondence between phase transitions in (tTc, 1) or 
[vTc, 1) and bad critical points p of ip, with the phase transition occurring 
at TT = \/ip{p). This includes vTc if and only if sup^g 11 V'(p) is attained and 
larger than A, in which case the last global maximum point is a bad critical 
point; if this supremum is finite not attained, then there is another first- 
order phase transition at tTc, while if the supremum is infinite, then tTc = 0. 
Finally, there may be a further phase transition at vr = 1 (with p = 1); this 
happens if and only if (/j(l) = /i(l) = A and 9?'(1) < 0. 

The phase transitions are first-order when the corresponding p is a bad 
local maximum point of (p, i.e., a bad critical point that is a local maximum 
point. (This includes tTc when sup(o,i] ^ is attained, but not otherwise.) 
Thus, the phase transition that occur are typically first order, but there are 
exceptions, see Examples 14.131 and 14.181 

Remark 4.8. The behaviour at vr = tTc depends on more detailed properties 
of the degree sequences (d^ )", or equivalently of Dn- Indeed, more precise 
results can be derived from Janson and Luczak [ill . Theorem 3.5], at least 
under somewhat stricter conditions on {di" )^] in particular, it then follows 
that the width of the threshold is of the order n~^'^, i.e., that there is a 
sequence tTcu depending on (dj )", with vTcn — > vTc, such that G*(n,d)7r,v 
and G*(n,d)^^e w.h.p. have a non-empty k-core if vr = VTcn + uj{n)n~^''^ 
with u){n) -^ 00, but w.h.p. an empty k-core if vr = Tr^n — w(n)n~^'^, 
while in the intermediate case vr = iTcn + cn~^'^ with — 00 < c < 00, 
P(G*(n, d)^ has a non-empty A;-core) converges to a limit (depending on c) 
in (0, 1). We leave the details to the reader. 

The same applies to further phase transitions that may occur. 

Remark 4.9. If A; = 2, then (jTS]) yields 

ip{p) := h{p)/p^ = ^Pjj(l - {f-py'^)/p, 

J>2 
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which is decreasing on (0, 1] (or constant, when ¥{D > 2) = 0), with 

sup ip{p) = hm ip{p) = y Pjjij — 1) = KD{D — 1) < cxd. 
P6(0,l] P^° j 

Hence 

vr^^^ = X/ED{D - 1) =ED/ED{D - 1), 

coinciding with the critical value in (13. lip for a giant component. 

Although there is no strict implication in any direction between "a giant 
component" and "a non-empty 2-core", in random graphs these seem to 
typically appear together (in the form of a large connected component of 
the 2-core), see Appendix |A] for branching process heuristics explaing this. 

Remark 4.10. We see again that the results for site and bond percolation 
are almost identical. In fact, they become the same if we measure the size 
of the k-core in relation to the size of the percolated graph G*{n, d),^, since 
v{G*{n, d)^^e) = n but v{G*{n, d).„y) ~ Bi(n, vr), so v{G*{n, d).„^y)/n -^ vr. 
Again, this is heuristically explained by the branching process approxima- 
tions; see Appendix lAl and note that random deletions of vertices or edges 
yield the same result in the branching process, assuming that we do not 
delete the root. 



Proof of Theorem. \4.5\ The case ¥{D > A;) = is trivial; in this case h{p) = 
for all p and vTc = so (i) applies. Further, Proposition 14.1( 1) applies to 
G*(n, d), and the result follows from Corefe(G'*(n, d),r,v) ^ Corefc(G*(n, d)). 
In the sequel we thus assume P(D > k) > 0, which implies h(p) > and 
hi{p) > for <p < 1. 

We apply Proposition 14. II to the exploded graph G*{h, d), recalling (j4.ip . 
For site percolation, G*(n, d),r,vi Pj = C~^'^Pj ^or j > 2 hy (j2.15p . and thus 

F{Dp=j)=C'7rF{Dp = j), j>2, 
and, because k>2, 

h{p) := E{Dpl[Dp > k]) = C'-^Hp), (4.10) 

hip) := ¥{Dp >k)= C^TThip). (4.11) 

Hence, the condition h{p) > Xp^ can, using (|2.18p . be written 

TTh{p) > a/. (4.12) 

If vr < TTc, then for every p £ (0, 1], by (fOj) . vr < tTc < Xp'^/h{p) so (|4.12p 
does not hold and h{p) < Xp^. Hence Proposition 14. l( i) applies to G*(n, d), 
which proves (i); note that if X^^Li ^ = 0{n) for some a > 0, then also 

h n n 

i=l i=l i=l j>l 

If vr > TTc, then there exists p G (0, 1] such that vr > Xp^ /h{p) and thus 
(I4.12P holds and Proposition I4.1( ii) applies to G*(n,d). Moreover, p in 
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Proposition 14. 1 ( ii) is the largest p < 1 such that ()4.12p holds. Since 7r/i(l) < 
^(1) < A and h is continuous, we have equality in (|4.12p for p = p, i.e., 
7rh{p) = XpF, so p is as asserted the largest p <1 with h{p)/{\p'^) = vr"^. 

Further, if p is a local maximum point of h{p) — \p^ = C,~^{-Kh{p) — Xp^), 
then Trh{p) — \p^ < in a neighbourhood oip and thus h{p)/{Xp^) < I/tt = 
h{p)/{Xp^) there; thus p is a local maximum point of h{p)/{Xp'^). Excluding 
such points, we obtain from Proposition 14.1( 11) using (j4.ip . (j2.14p . (j2.15p . 
dHHD, (|ilO]l and gUD, 



f (Core^) _ h t;(Core^) 
n n n 



(hi{p) = TThi{p), 

^cnDp=j) = 7^nDp=j), j>k, 



n n n 

e(Core^) h e(Core^) p Xp^ XpP' 

n n h 2 2 

This proves the result when p is not a local maximum point of h{p)/{Xp^). 

In particular, since hi{p) > 0, Core^ is non-empty w.h.p. when vr > vTc is 

not a local maximum value of h{p)/{Xp^). Finally, even if vr is such a local 

maximum value, we can find vr' with vTc < vr' < vr that is not, because by 

Remark 14.71 there is only a countable number of exceptional vr. By what we 

just have shown, G*(n,d)7r',v has w.h.p. a non-empty fc-core, and thus so 

has G*(n, d)^^v 5 G*(n, d)^/^v D 



Proof of Theorem \4.6\ We argue as in the proof just given of Theorem 14.61 
again using (j4.ip and applying Proposition l4.1l to the exploded graph G*{n, d) 
We may again assume P(L) > fc) > 0, and thus h{p) > and hi{p) > for 
< p < 1. We may further assume vr > 0. 

The main difference from the site percolation case is that for bond per- 
colation G*(n,d)^,e, ([2:27^ yields 

J > 2, 
J > 2, (4.13) 

h{p) := E(^pl[4 > k]) = (-'hipn'/''), (4.14) 

hi{p) := F{Dp >k) = C-^/ii(pvrV2). (4.15) 

Consequently, the condition h{p) > Xp^ can, using (|2.18p . be written as 
hip-K^f"^) > Xp"^, or 

TThipTT^/^) > A(pvri/2)2. ^^_^g^ 

If vr < VTc, then for every p G (0, 1] we have pn^'"^ S (0, 1] and thus by 

m 





F{D=j)=C-'F{D^,;,=j)., 


and hence 






F{Dp = j) = C'F{Dp^,,,=j 


and thus 





vr < vrc < 



/i(pvrV2) 
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SO (j4.16p does not hold and h{p) < Xp^. Hence Proposition 14. ll fi) applies to 
G*{h, d) as in the proof of Theorem 14. 51 

If vr > TTc, then there exists p G (0, 1] such that vr > Xp^ /h{p) and, as 
before, p is the largest such p and satisfies h{p)/{XjP) = vr"-^. Furthermore, 
if TT^/^ <p<l, then 

TTh{p) < TTh{l) < vrA < Ap^ (4.17) 

and thus p ^ P- Hence p < vr^'^. Let po := p/vr^''^. Then po G (0, 1] and 
Po is the largest p < 1 such that (j4.16p holds; i.e., the largest p < 1 such 
that h{p) > Xp^. We thus can apply Proposition I4.1l fii) to G*(n,d), with 
p replaced by po, noting that if po is a local maximum point of h(p) — Ap^, 
then p is a local maximum point of 

h(pn~y^) - A(p^-V2)2 ^ ^-l(/j(p) _ ,,-iAp2) 

and thus of 7r/i(p) — Ap^, which as in the proof of Theorem 14.51 implies that p 
is a local maximum point of h{p)/{Xp'^). (The careful reader may note that 
there is no problem with the special case po = 1, when we only consider a 
one-sided maximum at po: in this case p = vr^'^ and 7r/i(p) = Ap^ = Avr so 
h{p) = X and p = 1, vr = 1.) Consequently, when p is not a local maximum 
point of h{p)/{Xp^), Proposition SlKii) yields, using UO^Eh . KW . KT5\i . 

v{Covel) _ n v{CoTel) p 



n n n 



.cnDp,=j)=nDp=j), j>k, 

n n n 

ejCorel) _ n e{Coiel) p Ap| _ Ap^ 
n n h 2 2tt 

The proof is completed as before. D 

Consider now what Theorems 14.51 and 14.61 imply for the fc-core as vr in- 
creases from to 1. (We will be somewhat informal; the statements below 
should be interpreted as asymptotic as n ^ oo for fixed vr, but we for sim- 
plicity omit "w.h.p." and " — >".) 

If A; = 2, we have by Remark 14.91 a similar behaviour as for the giant 
component in Section [3l in the interesting case < tTc < 1, the 2-core 
is small, o(n), for vr < vTc and large, Q{n), for tt > tTc, with a relative size 
/ii(p(7r)) that is a continuous function of tt also at tTc and analytic everywhere 
else in (0, 1), cf. Theorem 13. Ill 

Assume now A: > 3. For the random graph G(n,p) with p = c/n, the 
classical result by Pittel, Spencer and Wormald [l?] shows that there is a 
first-order (=discontinuous) phase transition at some value c^; for c < Ck the 
/c-core is empty and for c > c^ it is non-empty and with a relative size ipk{c) 
that jumps to a positive value at c = Ck, and thereafter is analytic. We can 
see this as a percolation result, choosing a large A and regarding G(n, c/n) 
as obtained by bond percolation on G{n,X/n) with tt = c/X for c G [0,A]; 
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G{n,\/n) is not exactly a random graph of the type G*(n, d) studied in 
the present paper, but as said in the introduction, it can be treated by our 
methods by conditioning on the degree sequence, and it has the asymptotic 
degree distribution D ~ Po(A). In this case, see Example 14. 1 II and Figured! 
ip is unimodal, with ip{0) = 0, a maximum at some interior point po S (0, 1), 
and ip' < on {po, 1). This is a typical case; if has these properties for many 
other degree distributions too (and k > 3), and these properties of c/9 imply 
by Theorems 14.51 and 14.61 that there is, provided ^{po) > A, a first-order 
phase transition at tt = tTc = X/ipipo) where the fc-core suddenly is created 
with a positive fraction /ii(po) of all vertices, but no other phase transitions 
since /ii(p(vr)) is analytic on (tTc, 1). Equivalently, recalling Remark 14. 7^ we 
see that po is the only bad critical point of (p. 

However, there are other possibilities too; there may be several bad criti- 
cal points of (f, and thus several phase transitions of ip. There may even be 
an infinite number of them. We give some examples showing different pos- 
sibilities that may occur. (A similar example with several phase transitions 
for a related hypergraph process is given by Darling, Levin and Norris [d].) 




Figure 1. (p{p) = h{p)/p'^ for D ~ Po(lO) and k = 3. 



Example 4.11. A standard case is when D ~ Po(A) and k > 3. (This 
includes, as said above, the case G{n, X/n) by conditioning on the degree 
sequence, in which case we recover the result by [171].) Then Dp ~ Po(Ap_) 
and a simple calculation shows that h{p) = XpF(Po{Xp) > A: — 1), see [10|, 
p. 59]. Hence, if c^ := min^>o/i/P(Po(/i) > A: — 1) and A > c^, then 
TTc = info<p<i(Ap^//i(p)) = Cfc/A. Moreover, it is easily shown that h(jp)/p'^ 
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is unimodal, see [lO|, Lemma 7.2] and Figure [TJ Consequently, there is as 
discussed above a single first-order phase transition at vr = Cfc/A [171]. 

Example 4.12. Let fc = 3 and consider graphs with only two vertex degrees, 
3 and m, say, with m > 4. Then, cf. (j4.4p . 



/l(p)=3p3lP(^p 
= 3p3P^ 
Now, let P3 := 1- 



Dpl[Dp < 2] 

l)p2(l 



m[m 



D = 

-pY 



m 

-2 



)• 



3 \ D = 3) + PmE{Rp 

+ pm {mp - mp{l - pY''^ - 

a/m and pm '■= a/m, with a > fixed and m > a, and let 
m — > oo. Then, writing h = hm, hm{p) -^ 3p^ + ap for p G (0, 1] and thus 

hmip) 



(Pm{p) ■■-- 



pz. 



^oo{p) ■■-- 



a 
3p+-. 
P 



Since (f'^oi'^) = 3 — a, we see that if we choose a = 1, say, then (/'^(l) > 0. 
Furthermore, then 99oo(l/4) = | + 4 > C/9oo(l) = 4. Since also v?m(l) = 
3p3 — mpm = 3p3 — a ^ ip'^{l), it follows that if m is large enough, then 
(^^(1) > but (^m(l/4) > c^m(l). We fix such an m and note that ip = (pm 
is continuous on [0,1] with ip{0) = 0, because the sum in (j4.2p is finite with 
each term 0{p^). 

Let po be the global maximum point of 93 in [0,1]. (If not unique, take the 
largest value.) Then, by the properties just shown, Po ^ and pQ ^ 1, so 
Po G (0, 1); moreover, 1 is a local maximum point but not a global maximum 
point. Hence, tt^ = A/(^(po) is a first-order phase transition where the 3- 
core suddenly becomes non-empty and containing a positive fraction hi(po) 
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of all (remaining) vertices. There is another phase transition at vr = 1. We 
have C/9(l) = /i(l) = A, but since ip'{l) > 0, if ^i := supp < 1 : (p{p) > ^p{l), 
then pi < 1. Hence, as vr ^ 1, p{ir) y pi and h{p(7r)) y hi{pi) < 1. 
Consequently, the size of the 3-core jumps again at vr = 1. 

For an explicit example, numerical calculations (using Maple) show that 
we can take a = 1 and m = 12, or a = 1.9 and m = 6. 

Example 4.13. Let fc = 3 and let D he a. mixture of Poisson distributions: 
F{D = j) =Pj = Y, Qi P(Po(Ai) = j), j > 0, (4.18) 

i 

for some finite or infinite sequences (qi) and (Aj) with q^ > 0, J2i Qi = ^ and 
Aj > 0. In the case D ~ Po(A) we have, cf. Example 14.111 Dp ~ Po(Ap) and 
thus 

h{p) =EDp- F{Dp = 1) -2F{Dp = 2) = Xp - Xpe'^P - {Xpfe'^P 
= (Ap)V(Ap), 

where f{x) := (l — (1 + x)e~^) /x. Consequently, by linearity, for D given 
by dUl]), 

h{p)=Y,qiiXip)''f{Xip), (4.19) 

i 

and thus 

i 

As a specific example, take Aj = 2* and qi = X^ = 2~^*, i > 1, and add 
go = 1 — X^.j>i Qi and Aq = to make Y^qi = 1. Then 



H={p) = Y.f{2'p). (4.21) 



i=l 



Note that f{x) = 0{x) and /(x) = 0{x~^) for < x < oo. Hence, the sum 
in ()4.2ip converges uniformly on every compact interval [(^,1]; moreover, if 
we define 

oo 

V'(x) := E ■/■(2'^)' (4-22) 

then the sum converges uniformly on compact intervals of (0, cxd) and \'^{p) — 
'4>{p)\ = 0{p). Clearly, ■(/; is a multiplicatively periodic function on (0,cxd): 
'ip{2x) = ip{x); we compute the Fourier series of the periodic function ipC^^) 
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on M and find V'(2^) = X^nL-oo fp{n)e'^'^^^y with, using integration by parts, 
/o 



tp{n)= I 'il;{2y)e-^™y dy = / J^ /(2J'2S')e-2™?^d2/ 



j=-oo 



/•i /'OO 

V / /(2-'+^)e-2™ydy= / /(2^)e-2-i"ydy 
/(x)2;" 



j=-oo 

'■°° , rl-T- 

^-27rin/ln2_^ 



j;ln2 

1 
In 2 



1 f °° 

_L / (1 - (1 + x)e-")x-2™/ln2-2 ^^ 



^ ' j^e-^x-^'^^^/^^^-Mx 



In2(27rin/ln2 + l 

r(l-27rin/ln2) 

In 2 + 27rin 



Since these Fourier coefficients are non-zero, we see that V' is a non- 
constant continuous function on (0, oo) with niultiphcative period 2. Let 
a > be any point that is a global maximum of "0, let b £ (a, 2a) be a point 
that is not, and let Ij := [2~^~^b, 2~^b]. Then ip attains its global maximum 
at the interior point 2~^a in Ij, and since v'(p) — ''Pip) = 0{2^^) for p G Ij, it 
follows that if j is large enough, then ip(2~^a) > inax{ip{2^^^^b),ip{2~^b)). 
Hence, if the maximum of ip on Ij is attained at pj G Ij (choosing the largest 
maximum point if it is not unique), then, at least for large j, pj is in the 
interior of Ij, so pj is a local maximum point of (p. Further, as j — > oo, 

:p{pj) -^ maxip > V^(0) = l/ln2 while X:=ED = Y.i Qi\ = T.T '^'' = 1> 
so ip{pj) > A for large j. Moreover, Pj/2 G Ij+i, and since (j4.2ip implies 

oo cx> 

ip{p/2) = Y, fC^'-'p) = E ^(2V) > Vip), P>0, (4.23) 

i=l j=0 

thus (p{pj) < ip{pj/2) < ip{pj^i). It follows that if p G Ij for some i < j, then 
^{p) < ^iPi) < viPj)- Consequently, for large j at least, pj is a bad local 
maximum point, and thus there is a phase transition at ttj := X/cp^pj) G 
(0, 1). This shows that there is an infinite sequence of (first-order) phase 
transitions. 

Further, in this example ip is bounded (with sup(/3 = maxip), and thus 
TTc > 0. Since, by (j4.23p . supc/9 is not attained, this is an example where the 
phase transition at tt^ is continuous and not first-order; simple calculations 
show that as vr \ tTc, ^^(vr) = 0(7r — tTc) and /ii(p(7r)) = ©((tt — vTc)^). 

Because of the exponential decrease of |r(2:)| on the imaginary axis, |'0('T')| 
is very small for n 7^ 0; we have |0(±1)| ~ 0.78 • 10~^ and the others much 
smaller, so ^p{x) deviates from its mean ip{0) = l/ln2 ~ 1.44 by less than 
1.6 • 10~^. The oscillations of ip and if are thus very small and hard to 
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observe numerically or graphically unless a large precision is used. (Taking 
e.g. Xi = 10* yields larger oscillations.) 

Note also that in this example, (p is not continuous at p = 0; (p{p) is 
bounded but does not converge as p\0. Thus h is not analytic at p = 0. 

Example 4.14. Taking Aj = 2* as in Example 14.131 but modifying qi to 
2(^^2)1 £qj. gQj^g small e > 0, similar calculations show that p^tp{p) = '<p£{p) + 
0{p) for a non-constant function ^^ with multiplicative period 2, and it 
follows again that, at least if e is small enough, there is an infinite number 
of phase transitions. In this case, ip{p) — > oo as p — > 0, so tTc = 0. 

Since c/9 is analytic on (0,1], if there is an infinite number of bad crit- 
ical points, then we may order them (and 1, if (/^'(l) < and ip{i) = A) 
in a decreasing sequence pi > P2 > • • • , with pj -^ 0. It follows from 
the definition of bad critical points that then ^{pi) < ^{p2) < ■ ■ ■ , and 
supo<p<i v(p) = supj(/9(pj) = liuij^ooViPj)- Consequently, if there is an 
infinite number of phase transitions, they occur at {7rj}f° U {tTc} for some 
decreasing sequence vtj \ vTc > 0. 

Example 4.15. We can modify r4. 131 and 14. 141 and consider random graphs 
where all vertex degrees di are powers of 2; thus D has support on {2*}. 
If we choose P2' ~ 2~^* or P2i ~ 2^^~^'^ suitably, the existence of infinitely 
many phase transitions follows by calculations similar to the ones above. 
(But the details are a little more complicated, so we omit them.) A similar 
example concentrated on {10*} is shown in Figure [2l 

Example 4.16. We may modify Example l4.13l bv conditioning D on D < M 
for some large M. If we denote the corresponding h and (/? by Km and (pM, it 
is easily seen that Hm — > h unformly on [0,1] as M — > oo, and thus ipM -^ V 
uniformly on every interval [6, 1] . It follows that if we consider N bad local 
maximum points of ip, then there are N corresponding bad local maximum 
points of ipM for large M, and thus at least N phase transitions. This 
shows that we can have any finite number of phase transitions with degree 
sequences {dl"' )i where the degrees are uniformly bounded. (Example 14.171 
shows that we cannot have infinitely many phase transitions in this case.) 

Example 4.17. In Examples 14. 131 and 14. 141 with infinitely many phase tran- 
sitions, we have ED^ = X^i'?*!-^? + -^«) ~ ^^- This is not a coincidence; in 
fact, we can show that: // E D'^ < co, then the k-core has only finite number 
of phase transitions. 

This is trivial for k = 2, when there never is more than one phase transi- 
tion. Thus, assume A; > 3. If /c = 3, then, cf. ()4.4p . 

-2 



v{p) = h{p)/p' = Y, ^pi {p - p(i - pY^' - (' - i)/(i - py'^)p~ 

l>3 



/>3 
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Each term in the sum is non-negative and bounded by Ipiyl — {^—p))/p < 
lpi{l — 1), and as p ^ it converges to lpi{l — 1 — (^ — 1)) = 0. Hence, by 
dominated convergence, using the assumption YlviK^ ~ ^) — ^D{D — 1) < 
oo, we have ^{p) — > as p — > 0. For k > 3, h{p) is smaller than for k = 3 
(or possibly the same), so we have the same conclusion. Consequently, (p 
is continuous on [0,1] and has a global maximum point pQ in (0,1]. Every 
bad critical point has to belong to [pO)l]- Since if is analytic on [pojl]) it 
has only a finite number of critical points there (except in the trivial case 
<p{p) = 0), and thus there is only a finite number of phase transitions. 

Example 4.18. We give an example of a continuous (not first-order) phase 
transition, letting D he a mixture as in Example 14. 131 with two components 
and carefully chosen weights qi and 52 • 

Let / be as in Example 14.131 and note that /'(x) ~ 1/2 as x — > and 
f'{x) r^ — x~^ as X — > 00. Hence, for some o, A G (0, 00), ^ < /'(x) < 1 
for < X < 4a and ^x^^ < -f'[x) < Ix''^ for x> A. Let /i(x) := /(^x) 
and /2(x) := f{ax). Then, f[{x) < for x > 1. Further, if g{x) := 
/2(^)/l/l(^)l = {a/A)f'iax)/\f'iAx)\, then 

, , af'(a) a 



, . U K^Uj U, '-t 



af(4a) a/4 

A\n^A)\ > A.2(4^)- 

and thus ^(1) < ^(4). Further, if x > A/a, then /2(x) < and thus g{x) < 0. 
Consequently, sup^>]^g(x) = maxi<^<^/Q (7(x) < 00, and if xq is the point 
where the latter maximum is attained (choosing the largest value if the 
maximum is attained at several points), then 1 < xq < 00 and g{x) < g(xo) 
for X > Xq. Let /3 := ^(xo) and 

V'(x) := Pfiix) + /2(x) = pf{Ax) + f{ax). (4.24) 

Then ^'(x) < for x > 1, ^ijj^xq) = and %Ij'[xq) < for x > xq. 

Let b be large, to be chosen later, and let D be as in Example 14.131 with 
qi := /3aV(/3a2 + A^), q2 := I - gi, Ai := bA, A2 := ba. Then, by (lOOD . 

VP(p) = qi{bA)^f{bAp) + q^ibaffibap) = ±0^i;{bp). (4.25) 

Hence, ip'{xQ/b) = 0, ip'{x) < for x > 1/6 and ip'{x) < for x > xo/6. Con- 
sequently, xo/b is a critical point but not a local maximum point. Further- 
more, ip{xo/b) = qib'^A'^ f{AxQ)+q2b'^a'^ f{axo) and A := ED = qiXi+q2^2 = 
b{qiA + q2a); hence, if b is large enough, then (/?(xo/6) > A. We choose b such 
that this holds and b > xq; then p := xo/b is a bad critical point which is an 
inflection point and not a local maximum point. Hence there is a continuous 
phase transition at fr := X/ip{p) £ (tTc, 1). 

We have ^'{p) = ^"{p) = 0; we claim that, at least if A is chosen large 
enough, then (p"'{p) / 0. This implies that, for some 01,02,03 > 0, ip{p) — 
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(p{p) ~ —Ci{p — p)^ as p —>■ p, and ^(vr) — p(vr) ~ C2(vr — tt)^'^ and hi{p{TT)) — 
hi{p{TT)) ~ cs{tt — tt)^'^ as vr ^ vr, so the critical exponent at fr is 1/3. 

To verify the claim, note that if also ip"'{p) = 0, then by (|4.25p and (|4.24p . 
'(/''(a^o) = '^"{xq) = '4j"'{xq) = 0„ and thus 

/3AV(^')(>lxo) + aif^^\axQ) = 0, j = 1,2,3. (4.26) 

Let xi := ^xq and X2 '■= clxq. Then xi > A, and f'{x2) > so X2 < C for 
some C. Further, (j4.26p yields 

X2f"{x2) _ Xir(xi) Xi/-"(X2) _ Xfr(xi) 

/'(X2) /'(Xl) /'(X2) /'(Xl) • ^ • ^ 

Recall that xi and X2 depend on our choices of a and A, and that we always 
can decrease a and increase A. Keep a fixed and let A — > oo (along some 
sequence). Then xi — > oo but X2 = 0(1), so by selecting a subsequence we 
may assume X2 — > y > 0. As x — > co, /'(x) ~ — x~^, /"(x) ~ 2x~^, and 
f"'{x) ~ — 6x~^. Hence, if ()4.27p holds for all large A (or just a sequence 
A — > oo), we obtain by taking the limit 

^^^^ = hm :^:^^ = -2 and ^^^^^ = lim ^^^^^ = 6. 
f'iy) ^^°° fix) f'{y) a^^oo /'(x) 

Finally, let F{x) := x/(x) = 1 - (1 + x)e-^. Then F"{y) = yf"{y) + 
2f{y) = and F"'{y) = yf'jy) + 3/"(y) = 6y-'ny) - 6y~'f{y) = 0. On 
the other hand, F'{x) = xe~^, F"{x) = (1 — x)e~^, F"'{x) = (x — 2)e~^, 
so there is no solution to F"{y) = yF"'{y) = 0. This contradiction finally 
proves that (p"'{p) 7^ 0, at least for large A. 

5. Bootstrap percolation in random regular graphs 

Bootstrap percolation on a graph G is a process that can be regarded as a 
model for the spread of an infection. We start by infecting a subset ^0 of the 
vertices; typically we let ^0 be a random subset of the vertex set V{G) such 
that each vertex is infected with some given probability g, independently of 
all other vertices, but other choices are possible, including a deterministic 
choice of ^o- Then, for some given threshold £ E N, every uninfected vertex 
that has at least I infected neighbours becomes infected. (Infected vertices 
stay infected; they never recover.) This is repeated until there are no further 

infections. We let Aj = Ao be the final set of infected vertices. (This is 
perhaps not a good model for infectious diseases, but may be reasonable as 
a model for the spread of rumors or beliefs: you are skeptical the first time 
you hear something but get convinced the £th time.) 

Bootstrap percolation is more or less the opposite to taking the /c-core. 
For regular graphs, there is an exact correspondence: it is easily seen that 
if the common vertex degree in G is d, then the set V{G) \ AK of finally 
uninfected vertices equals the (d + 1 — £)-core of the set V{G) \ Aq of ini- 
tially uninfected vertices. Furthermore, if the initial infection is random, 
with vertices infected independently with a common probability q, then the 
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initial infection can be seen as a site percolation, where each vertex remains 
uninfected with probability n = I — q. Consequently, in this case we ob- 
tain results on the size of the final uninfected set from Theorem 14. 5^ taking 
k = d — i + 1 and tt = I — q. 

Bootstrap percolation on the random regular graph G{n,d) with fixed 
vertex degree d was studied by Balogh and Pittel [l| . We can recover a large 
part of their results from Theorem 14. 51 We have, as just said, k = d — i + 1 

(n) 

and TT = 1 — q. Moreover, all degrees d) = d; hence the definitions in 
Section [2] yield rij = ndjd, Pj = 5jd, £>„ = d, D = d and X = KD = d. 
Condition 12.11 is satisfied trivially. Furthermore, (|4.2p and ()4.3p yield, since 
7^p~Bi(d,p), 

Hp) = JZ^^dM = JZ^['^)p^{i-p)''~' = E^^f^^I l)^"'(i -^)'"' 
j=k j=k ^-^^ j=k ^^ ^ 

= dp¥{Bi{d - l,p) >k-l)= dpF{Bi{d -l,p)>d-i) 

= dp¥{Bi{d-l,l-p) <i-l) 
and 

hi{p) =F{B[{d,p) >k) =P(Bi(d,p) >d-£ + l) =P(Bi(d,l-p) <£-l). 
Consequently, ()4.9p yields 



TTr := mi . , , = mt 



o<p<i h{p) o<p<i F{Bi{d -l,l-p) <£-!)' 

We define qc := 1 — tTc, and Theorem 14.51 translates as follows. (Recall that 
we have proven Theorem 14.51 for the random multigraph G*{n, d), but as said 
in the introduction, the result holds for the simple random graph G{n, d) by 
a standard conditioning.) 



Theorem 5.1 (|l|]). Let d, £ and q £ [0,1] be given with 1 < £ < d — 
1. Consider bootstrap percolation on the random d-regular graph G{n,d), 
with threshold £ and vertices initially infected randomly with probability q, 
independently of each other. Let 

n = Q^^^ •= 1 — inf - ib 1) 

■ o<p<iP(Bi(d-l,l-p) <£-l)' ^'' 

(i) If q > qc, then \Af\ = n — Op{n). Furthermore, if I < d — 2 then 
w.h.p. \Af\ = n, i.e., all vertices eventually become infected. 

(ii) If q < qc, then w.h.p. a positive proportion of the vertices remain 
uninfected. More precisely, if p = p{q) is the largest p < 1 such that 
P(Bi(d-l,l-p) <£-l)/p = {l-q)-^, then 

|^/|/n^l-(l-g)P(Bi((i,l-p) <£-l) < 1. 

Proof. It remains only to show that in case (ii) , p is not a local maximum 
point of (p{p) := h{p)/{dp^) = P(Bi(d- 1, l-p) < £-l)/p. (In the notation 
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of Section HI (p{p) = ip{p)/X.) In fact, some simple calculus shows, see [l|, 
§3.2, where R{y) = (p{y)~^] for details, that the function ip is unimodal when 
i < d — 1 and decreasing when i = d — 1; thus there is no local maximum 
point when i = d — 1, and otherwise the only local maximum point is the 
global maximum point po with (^(po) = ^^c^ = (1 ~ 'lc)~^- (It follows also 
[]J that the equation (p{p) = {1 — q)^^ in (ii) has exactly two roots for every 
q < Qc when I < d — 1 and one when i = d — 1.) D 

Remark 5.2. The case i = 1 {k = d) is rather trivial; in this case, Af 
is the union of all components of G(n, d) that contain at least one initially 
infected vertex. If further d > 3, then G{n, d) is w.h.p. connected, and thus 
any non-empty ^o w.h.p. yields \Af\ = n. (The case d = 2 is different but 
also simple: G(n, 2) consists of disjoint cycles, and only a few small cycles 
will remain uninfected.) 

Actually, Balogh and Pittel []| study primarily the case when the initially 
infected set Aq is deterministic; they then derive the result above for a 
random Aq by conditioning on Aq. Thus, assume now that Aq is given, 
with l^ol = '>^- (For G{n,d), because of the symmetry, it does not matter 
whether we remove a specified set of m vertices or a uniformly distributed 
random set with m vertices.) Assuming m ~ nq, we have the same results 
in this case, see Remark ll.lt indeed, the proof is slightly simpler since the 
use of the law of large numbers in Subsection 12.11 is replaced by the obvious 
hd = n — m, hi = rij^ = dm. 

Theorem 5.3. Theorem \5. 1\ remains valid if the initially infected set is any 
given set with m = m(n) vertices, where m/n -^ q. 

Remark 5.4. As in Remark 14.81 it is also possible to study the threshold 
in greater detail by allowing q to depend on n. If we assume i < d — 2 (i.e., 
A; > 3) and q = q{n) —>■ q^ defined by (j5.ip . then Janson and Luczak [ill . 
Theorem 3.5] applies and implies the following, also proved by Balogh and 
Pittel jl| by different methods. 

Theorem 5.5 ([ifl). Consider bootstrap percolation on G{n, d) with i < d—2, 
and assume that the set Aq of initially infected vertices either is deterministic 
with \Ao\ = q{n)n or random, with each vertex infected with probability q{n). 

(i) If q{n) — ^c S> n^^i"^ , then w.h.p. \Af\ = n, i.e., all vertices become 

infected. 
(ii) If qc — q{n) ^ n^^''^ , then w.h.p. \Af\ < n and, moreover, 

\Af\=hi{p{q{n)))n + Op{n'/^\q{n)-q,\-'/^). 

Janson and Luczak [11, . Theorem 3.5] is stated for random multigraphs, 
and for G*{n,d) it yields further an asymptotic normal distribution in case 
(ii), as well as a precise result for P(|^j| = n) in the transition window 
q{n) — qc = 0(n~^'^). The latter result can easily be transformed into the 
following analogue of [ill . Theorem 1.4]; the asymptotic variance o"^ is given 



by explicit but rather complicated formulas in 111 ]. 
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Theorem 5.6. Assume i < d — 2. Infect (from the outside) the vertices 
in G* (n, d) one by one in random order, letting the infection spread as 
above to every vertex having at least i infected neighbours, and let M be 
the number of externally infected vertices required to make \Af\ = n. Then 

(M - nq^jn^l'^ -^ A^(0, o'^), with a^ > 0. 

Presumably, the same results hold for G{n,d), but technical difficulties 
have so far prevented a proof, cf. [llj. In any case, it follows from Theo- 
rem [5]6] that the size of the transition window is 0{n^^''^) for G{n,d) too, 
and not smaller. 

Appendix A. The /c-core and branching processes 

We give a precise statement of the relation between Proposition 14.11 and 
branching processes. This can be seen heuristically from the branching 
process approximation of the local exploration process, but as said above, 
we do not attempt to make this approximation rigorous; instead we compare 
the quantities in Proposition 14.11 with branching process probabilities. 

Theorem A.l. Let X be a Galton-Watson branching process with offspring 
distribution D* and starting with one individual o, and let X be the modified 
branching process where the root o has offspring distribution D but everyone 
else has offspring distribution D* . We regard these branching processes as 
(possibly infinite) trees with root o. Further, let T^ be the infinite rooted tree 
where each node has k — 1 children, and letT}^ be the infinite rooted k-regular 
tree where the root has k children and everyone else k — 1. 

Then p = T{X 5 7^), the probability that X contains a rooted copy of Tj. 
(i.e., a copy of T^ with root o) and hi{p) = W{X D T^), the probability that 
X contains a rooted copy ofTj.. 

Hence, by Proposition 14.11 the probability that a random vertex be- 
longs to the fc-core, which is E(f (Core^)/n), converges to the probability 
P(A' D Tfc), the probability that the branching process approximating the 
local structure at a random vertex contains the infinite fc-regular tree T^. 
Similarly, the probability that a random edge belongs to the Core^, which 
is ~ E(e(Core^)/(nA/2)), converges to h{p)/\ = p^ = ^{X D Tkf, which 
can be interpreted as the probability that both endpoints of a random edge 
grow infinite A;-regular trees in the branching process approximation. 

Proof. Let T^n be the subtree of 7^ consisting of all nodes of height < n, 
i.e., the rooted tree of height n where each node of height < n has k — 1 
children, and let qn be the probability that X contains a copy of 7fc„- Thus, 
go = Ij and for n > 0, (?n+i is the probability that the root o has at least 
k — 1 children that each is the root of a copy of T^n in the corresponding 
subtree of X; let us call such children good. By the branching property, the 
subtrees rooted at the children of o are independent copies of X: thus the 
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probability that a given child is good is g„, and the number of good children 
of o has the tinned distribution D* . Hence, 

oo oo 

(?„+i = ¥{Dl >k-l) = Y, nD* =d-l)Y^ JP(Bi(^ -l,qn)=l- l) 

d=k l=k 

d>k l>k ^" 

=i|:'r(o,.=o=^'.(,.). 

Since x i-^ h{x)/{Xx) is increasing (e.g. by the same calculation) and 1 = 
Q'o ^ 9i ^ • • • ) it follows that qn decreases to the largest root p of j^h{q) = q 
in [0,1]. On the other hand, the events £n := {^ 5 '^kn} are decreasing, 
^1 2 '^2 5 • • • ) and Pl^ £n is, by a compactness argument, equal to the event 
{X 2 %}■ Hence, F{X D %) = lim„ q^ =_p. 

Similarly, X contains a rooted copy of T^ if and only if the root o has at 
least k good (now with n = co) children. We have shown that each child 
is good with probability p, and thus the number of good children has the 
thinned distribution Dp; hence F(X D T^) = F{Dp > k) = hi{p). D 

Remark A. 2. When k = 2, T2 is just an infinite path, and thus p = 
¥{X 5 T2) is just the survival probability p of the branching process X, as 
observed algebraically in Remark 14. 4[ Hence the thresholds for 2-core and 
giant component coincide, for any of our percolation models. Moreover, we 
see that if t> is a random vertex, the events "u is in a giant component" and 
"u is in the 2-core" are approximated by "the root o in ^ has at least one 
child with infinite progeny" and "the root o in X has at least two children 
with infinite progeny" , respectively, which again shows the close connection 
between these properties. 
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